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۶ قاطا ال 0 ۲۱۵۲ 11-11-01 :۹۵661۷60 
0ج مها تصحطحلوط دول ریق تهمطه ۷۵ بیط طقصمم‌صه ۷22 یو ۲۵8۵221 13-12-01 :1۱6۷1660 
عصلون] وتمط0۳0 م۸ 10۲ عمتصنه ۷۷ «اتقط ۳۲۵۹۲ عصزیمد .2022 ۸۲۰ ۱۷۲۵۲۵۸012۵۸06 28-02-02 :60660۲60 
اوعد طز۳ جهتوزع۳ م1) .127-144 :(36)01 9011 هه ۲ع2 ۷۷ ۵۶ لقصنتامز .۱۵061 ۲ ۱۷۷ 20-05-02 :0۵1106 ۷۵11201 
(2050201 
7 ۷ :]10[ 


م۱۹۵۵ 

۴۲ 1 اهتدم مستتالته‌تتعة م6 وعقوهاً أَمتمصقصت پ«حقعط فعوممرصصا هصتیصو صا وم فامالميم احتتاانه‌تعخظ 
۷ 2 18 ]از 2150 اه ومامصصهتوم ماقصصتاه اصه)زمحرص اومجطر فص ۵۶ مهم و1 متااجهم‌ها چاصم )ولا .عل‌تقطهته و خصوت1 ۵۶ 
۷ 276 طمنطه ۶2۱1 هه و۲ مه فصمتامنهع تمصع «امصصتاصنا .ص۹۵6 آمتمانمتهه مفطا صا غصمصصماه متام 
ق0 تن لته معط ممتاحعلنفجمع ماصا عصت‌له) معط 1 .ومعمصصهل 62۷۷ 20 عممطاو مفییه ۱۷۱۱1 ۱۵۵2۵66060 
0 ممزلن0د م1 .دمم )۲۵۱۵6۷2 عمتمنتع؟ 10۲ ممصه)۲ممرح کهمتع ۵ و0 ۷۵ طمهمتومه مه عصتنع]1ه 2240 ]1۲05 
۶ وم مهن هه ممتافتم‌امهتهطل (2015) .2 6۲ زیهدع۲۱ مه (2012) 200۱0مظ مقطمطهنز مه عدبلتمه ها غنام 
۲ 10۵ 1800010660 ۷۷۵۲۵ 0615مجظ اممتم)ن۱ ,ولمم احعتانام‌تبوه معط ما ممتقاهع ما 06۲60تفصمن صعه 2۷۵ظ وو۲۵] 
۲ ۵06 ]و۵ عطا عصتصتصمام ما ما20 ما رومتلته محر مر .فصمنعع1 0ماهعناوع۷ 1۳ اصمتعنل طا مق فعطلنه10۳] 
۰ مه 60هصتاوه جمع‌ها فقط )وم افقا ۶ مق ماطوطام0۳ ,ععفاه عصته1۱۵۳۷ ۵۶ معصجسهعه ۵۶ ملق مطا عصتاه‌ن2۳60 
6۵ 220 (2017) ۲مجمونای [-2مص) نها رها 4عنامحرج حمتطه 0مطاعصه 2 و1 دععصقطه معبمهتهم‌جصعا حصعا همم مصلاهعنا۵ع 18۷ 
تفحآ .0216 ]1۳05 )فقاً عطا طاً مععصقطه مفلج مه مق عصتیه10۳ ما دععصقطه 0مطفناطاهاوه غامهاه عیام فص م (2018) .2 ام 
۸ 06 موه عممم۲۵ا06۵ 1۲۵9۲ 1۵۲ عصتصنه ز«اتقه سامط- 48 هط ۲۵211260 (2019) طقصهم‌مهل۷22 هه تصحطم]و۳۴ 
1 ,ما02 م۲۱۵۳ متس اتملنمتانهن ]وم منم ۵ 1610 مطا صا فمتنه 2 مالموعنا ,صمتفتمعتن متقومم20 طاز۳ 
.0هط0زه طمهه عم همه اعد معلقصه مه 0۷1060عج مضه 60حول1صاهاوع ۵ اوناظط مصمادره عصتصته۷ ]1۳05 ۲۵۵10 2 966125 
۷۵20 10 2۷6 ۷۷۶ 30 رعصتا1۳05 ۵۵10۲۵ درجل معط مصتصه اتمه طمیده معلقصظط مه و لمع عنام ممصله رجحقوع فنطا 
۵ ,200 تعطان فطع و اع0مصص ۱۷۷ تون ممتاهایصته معبهتمم‌جصهعا صیصصتصتمطر تتامج- 72 ۵۶ اتلتله۷ مه رمح‌تنامع2 
۵ 2960 و2227 ]و۳0 اباوماج ععصریع! فص مصررمصا ما مين طمقع صا امنلمتج طمقع 0۴ معهاه اهمزع0]مممصام 120۷ ]5ناظظ 
27 و ا۲0ظ) وععهاه لهمزعم1معمطم )06060 مه و1 62و۲6 فنطا گم لمع ۹66080 مط) ۵۲۵۶0۲ رمتنتاههم‌صما آهعتاتره 

220 1280526 8 ۰ 


۵۹ 2۱0 واوز۱۷]۸)۵ 

۵۵۲ ۶() دفباهعطا ممتح‌ااصته متنقعهم‌مها میمصت تیامط-772 رتاو اصمتستن ۵۶ صصتح عمط عمتطعج مه تع۵۳0 1 

0 2 10۲ 06۲1۷۵0 ۷۷۵۲۵ 1۳067 ط۷۵۵۵)۵0 ۵۶۴ ومل۷2 مه 18۷65022660 ۷۷25 2۵061 (ظ۷۷/۳۴) عصتا۳۲۵:۵۵28 204 ۵562۲60 
040 [2همناه] نامهم .2016-2107 هط ,۲۶۵۷186 صوزنه۸2۵:0 ۷۷۵۹۲ ,م2اتطهک ما هتمطمتن آه)صمجصتنهمنه مه 26 آملم وتعامظ 
۶ ممای امومع ۱۵۲209 ط مزع 0ع)ومظ معط ۵ فصتهاجم املممه ۷۷۴۲ عمتمیا ممتاهاناصصته ماجهم‌طها فتمام9 2 10۲ 
1 2000) 5660080 30 مه لقاع و1 0242 6وبا ص12 لمح )طونم صتحععا ۵۶ ممتاناموع: ۲۱۵۲۱20۵1۵ .وم)مطملن! 3 20 9 ,27 
مدرگ افدعع۲۵۲ لجطامات فطا حصم ممتاماموعد لمیم۵۲2۵ظ 0.5 ط1 فصمتاتهجمی تجل‌صنمه ط-3 0صه لهتانصا م1 .رح 
0۳6۷1۵۱۵ ما ج۵ 0عفظ ,رطلن۱) ممتامصصملم اممممصممطظ رم فاصم افممتاهآظ صم 0مصتفاهاه مهس (د0)) 
۳۱8062466 ۷۱۵0615 .ظمتهعوعع ونط مد ۵۵ونا مه [06مصه ۷۷ 10۲ ممتاماع‌تومن مصصمدامه لهم1وررطاص [ 1۳۱۷۲۷ جاهتهعوع۲ 
0 معمتامانم‌صم تقصماصاً صاصق ۵0 )هارمه گم فممتاهای تمطاهع ۳۵۵و 51 ۵۶ مه 10۲ تام 171 03:00 
فعنله۷ 0عصنهاطاه عطا صقط رل‌مطامصه تفمصنا مصنوها ممتاعصنگ 3 تماما مهو شا ۸ ها 0عاقم01معاظ ۷۸۵۲۵ 


,0602۵00 [2م1ووظظ ۵۶ )صمحصتدمون1 رو۶۵]6990۲ظ ملهلهموفض مه بعمامامصتات دوه ۵۶ اصعل‌تاو .۲۲۰ظظ -3 4صه 2 ,1 
۰ ,1۲۳۵1 رصهطه۶ع1 بممطهلوز ۵۶ «انونمنهرنا رعصتصصماظ مه فممممتمی مومع ۵۶ وانامج۲ 

(15. 8602.01.26 ۵) هل ۰۷۵2 :اتممصظ ممطه مصت‌طمصوعتزمن -۶) 

صحی رلعمامطفطه رالوهت نا 0عممطمطه ,اه‌هد له ولنمع۲ بممتلهعن۱ گه اصمانهم۱(6 ررمووع]۵:ظ تصماوزووظ -4 

۶ اهنا هعموص ممتاماامموصه:] مه 0 گه امد رومتامههمع0) ۵۶ اصمحصنهم 1(۵۵‏ ر990۲ع۲۳۵1 اصداولوع۸ -5 
0 ,1512۳2۳ مصفطع]15 


۸ نشریه آب و خاک. جلد ۰۳۶ شماره ۱. فروردین - اردیبهشت ۱۴۰۱ 


۷۵ 24 ۴۱۷۲9۲ ,۷۵ رها عصنوت رها فعمتاعاه مطا ما ۲۵0موهاه تمصع حصتتمتصتحظ مه 0عتقم‌طرهن ۰ ۱۷۸۵۲۵ 
صحیا حصم 0مصنماهاه ۷۵۲۵ وععهاه طاسمتع گه عم م4صه عصتصمتععه که مصصتا مطا روعتافتاهاه لمعنع0[معمطظ .1۳01620018 
5 8 )1052 41 4ص2 رتمقصمو +۳۲۲۳۷/۲ ۵۶ ومهمصطا منامام 7 1۵20521 77 ,5ع910عظ .صمتام2لصموع ل2ع1ع۱۷]۵)6۵۲010 
-2007 )560 ما طمته حصرم ماتوصام هبو لهمنوم۱هع0 ومیهاو: 0ع)تصنا حصم 0007۵1۵۵060 ۷/۵۲۵ ۲مقصعو نام 0۶ 
۵ 0۲۲۵۵10 . 10تامصطم۲201 مه متتمطامدمصصاه رطل‌تهعوع۴ فتطا ص .عتعامصه 30 که ممتانااووع۴ لفمتادمو 2 ۳16 2016 
معط مجح )ممصصصمتزصه متهسمو ۳۵۲۷۱ فطا طصا علق. م)ممامصه فطا عم 0مط)عصهر صهمضصاط فطه ماه مهم 

۱ ۱۱۱۷ عطروتا 6210012060 قه 1406 ممتلمامهعع۷6 


مکی وز(۱ 20 وااباععع 
40 ممتاهاهمه اصمه‌تروزه ۵ ۲۵۷۵۵۱۵0 فالناوع۲ راملمصر ۷۷۲۲ فط ۵۶ وعمهع401ص1 ممتاهتويه مطا عصتطتی؛۲ 
مج معا اتمه مرج ناوتان مضه ۱۷/۴ رم مامافتنه۷ ممتلجوم‌مصعا مماممط 2 06007۵۵ 120061 ۲۵۵1۲۵8910 
-72 10 2001160 06 صق م۵0 ۷۷۲۲۴۳۲ ااناوعع ه مه .صمتاقاناصصته عتامط-72 10۲ محمتاهاه معقاجه مط صا ۵0رموهان ما0ه۷۵2۲1 


0 طمولتهمعهه صاً خقط) 10162060 جم‌نهعوع: ونطا گه مصتهصظ تعطامصن .تاه ۵۶ هه معط ط ممتاها نامه متاجتهم‌جطها عتامظ 
0 وعوصفل ممهعم1۳4 زآهممتن ومعفمصظا اقفل‌صما مق عمتجم فمبتله۷ ۱۲۳۷ بع26 11010-76607060 عطا 
۵061 500۷7۵0 زمره 1390061 مص مامالهبه هه 0مویا قم6هع01ض1 مط) رتتاه ونط م1 تقطمعن ماومه )صونعام مطا طا دعهماه 
-72 عطا )اه مطمتاهاناصصته عتامح-72 ۲۵۲ صقطع عممتاهااصصته تام 48 صفطع 0صع عتام24 0۲ عهتناهعع عتمصهر منج 2665ع0طنظ 
+2 طع)8ع/۱ظ):۵ ما .عصمتاهانامصته عتامط 48 مضه عبامط-24 حصمی کصمععل طمننج عمط ود پامجتتاممه صمتهاتاحصته تام 
000۵۵0 مهن از فهعتة صا ب«ماتج مفقلندته م6 ]من م۷ و1 1 رجمتعع فبامصتقمتممه مه ود ت۷۲ 
۶۲ 36 20 ,0.5 ۶ معص۲2 مطا طاً و1 فصمتاهاصته معط الم ما مصمتماه له ۶ )همتم‌تگلههه ممتاجاعته لهاما مطا 1016ع7عظ) 
عنطا ۵۶ عصنص 96609۵ مط) ما مصللع0عع۸ .عباتواعن 3.8 20 2.8 )امه فعطمهع۲ ۲۵۹۵6۵۷۵ ,۴۳۲9 مه له ۵۶ و۲26 
0ص ووعتعمتم مطا «مطو صقن وععمحصصا ماتالمامه ۵ )82صض1 مضه 7 52صضنا حصم 0مصتقاهاه 10162001 ۲ ۱۲۱۷ عطا رطاهه6وع۲ 

۰ 2001۱6 0۴ ومعهاه 09۵01021621 مطاً مد وععصفطه 


جمتعصلصه) 
کج متااهزهم‌ججرما ممریتحتصتصر عص ۵۶ ممتتهلنصته عتامط-72 ۲۵۲ زملممه ۷۷۲ عمط هم پ«مممته‌تاه معط 0ع7مطو رنه وتط ]1 
1 007و عطا ما قمهماه لممنع10مصعطم مامرمه عصتامع۱ع0 طاً مععقصصا 8 )92صم مه 7 0921صضا ۵۶ لهتاجماه معط مه ۷۸۵۱1 
مط2نع۵عع مه ممتهایصته مها تور 72-0 ۲۵۲ آع0مصهر ۷۷ فطا. عصتفنا. ها رعت0]ع7عط][ 
2 ۲ 19 ما ۱۵۷۵ م0 2 وعطامهع 0عقطهته هچ صا معتتاجهموصصعا مطا 1۴ روععمصظا 052صقا ما ممعهاه [102162مصمتام 
۷ 0006 0۵6 10امطو طمتاهع‌تاتهه ]05 معط 4ج 0ععصیامصصج ع۵ صقع عصتصنه/۳ ]1۳05 رمعفاه لهمنع10مصمطام عطع م6 مب وعتامط 
۰ 38 


۵061 ۷۷۴۹۳ ,)۲05 عصترمگ رز ۷ ۱۱ روه‌عهاه لهمنعم1مصمطم گه ممنام ۲۱06۵ رماوره مه اه ۱۵۲۱۷۵۲ 


8 نشریه آب و خاک 


۷ 


جلد ۰۳۶ شماره ۱ فروردین- اردببهشت ۱۴۰۱ ص ۱۴۴-۱۲۷ 


پیش آگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۲۲1۴۲ 


سیده محبو به اين ححازی ۱- ححت اله یزدان پناه "- سعید موحدی۲- محمدعلی نصر اصفهانی *- مینا مرادی زاده؟ 


تاریخ دریافت: ۱۴۰۰/۰۸/۲۰ 


تاریخ پذیرش: ۱۴۰۰/۱۲/۰۹ 


چکیده 


سرمازدگی محصولات کشاورزی در فصل بهار همه‌ساله زیان‌های سنگین مالی را به بخش کشاورزی به‌ویژه در باغات شمال غرب ایران وارد 
می‌کند. هدف این مقاله, ارزیابی سیستمی برای پی شآگاهی سرمازدگی با استفاده از شبیه‌سازی دمای حداقل ۷۲ ساعته به‌وسیله مدل ۷۷1۳ و تشخیص 


مراحل فنولوژی سیب از تصاویر لندست است تا با شناخت مراحل فنولوژی محصول و دمای بحرانی در آن مرحله چنانچه دمای حداقل در ۷۲ ساعت 
آینده به دمای بحرانی برسد پیش‌آگاهی سرمازدگی صورت گیرد. داده‌های دمای ۲ متری خروجی مدل ۷۷18۳ برای شبکه محاسباتی داخلی» در ۵۱ 
ایستگاه سینوپتیک با دمای حداقل مشاهداتی در ایستگاه‌ها مقایسه شد. مقادیر شاخص ]2۷ نیز با استفاده از تصاویر لندست ۷ و ۸ سنجنده‌های 
+۲۲۱۷ و 1 در سال‌های ۲۰۰۷-۲۰۱۶ برای باغ سیب واقع در ایستگاه تحقیقات هواشناسی کشاورزی کهریز ارومیه محاسبه و با زمان مراحل 
فنولوژی ثبت‌شده در محل مقایسه شد. نتایج نشان داد که معنی‌داری همبستگی و مدل رگرسیونی بین متغیر دمای ۲ متری خروجی مدل ۷۷1۳ و 
متغیر دمای حداقل مشاهداتی در مجموع کل ایستگاه‌ها برای شبیه‌سازی ۷۲ ساعته وجود دارد. درنتیجه می‌توان از مدل ۷۷16 در شبیه‌سازی ۷۲ ساعته 


دما در منطقه موردمطالعه بهره برد. يافته دیگر این تحقیق نشان داد که در مقایسه با داده‌های زمینی ثبت‌شده در منطقه» مقادیر ۱12۷ به‌دست‌آمده از 
تصاویر لندست به‌خوبی گویای تغییرات مراحل فنولوژی در باغ سیب موردمطالعه است. 


واژه‌های کلیدی: تشخیص مراحل فنولوژیکی, سرمازدگی بهاره سیستم هشدار سریع. مدل ۱۷181 ۷9۷1 


مقدمه 

دما علاوه بر آن که یکی از پارامترهای مهم قلیمی است عنصری 
بسیار تآثیرگذار در بخش کشاورزی نیز می‌باشد. نوسانات و 
افت‌وخیزهای نابهنگام دمایی که اغلب دور از انتظار است باعث 
غافلگیری و بروز خسارت سنگین می‌شود ( ,.۵1 2 160912۷1 
8 هرگاه در فصل بهار که زمان رویش برگ‌های جدید و 


۸ ۲ و ۳- به‌ترتیب دانشجوی دکتری آب و هواشناسی کشاورزی و دانشیاران گروه 
جغرافیای طبیعی, دانشکده علوم جغرافیایی و برنامه‌ریزی» دانشگاه اصفهان» 
اصفهان» ایران 
(*- نویسنده مسئول: 
۴- استادیاره گروه مهندسی آب دانشکده کشاورزی, دانشگاه شهرکرد شهرکرد 
ات 
۵- استادیار. گروه مهندسی نقشه‌برداری» دانشکده مهندسی عمران و حمل و نقل؛ 
شاه اما استهات یاه 

10]: ۷ 7 


7۴ ۵ 0120 ۰۷۵2 :اتفح۲) 


پرطراوت است دما به زیر نقطه انجماد برسد خسارت سرمازدگی بهاره 
رخ می‌دهد (1993 ,اعقلم2 0صد ملناعاه:ذ۳1). در بیشتر مواقع این 
روی آسمان بدون ابر یا به دلیل سرد شدن لابه هوای نزدیک به 
سطح زمین درنتیجه تابیدن گرما از سطح زمین به جو اتفاقق می‌افتد. 
ماهیت خسارت سرمازدگی به مرحله‌ای از رشد گیاه بستگی دارد که 
سرمازدگی در آن رخ می‌دهد. به‌ طورکلی در ابتدای رشد و نمو گیاه 
مقاومت آن به سرمازدگی بسیار کم است ( ,۹260 4ص 0:22 
0 یکی از مراحل فنولوژیکی بحرانی در درختان میوه خزان 
کننده مرحله شکوفه دهی است زیرا در این مرحله گیاه به تنش‌های 
محیطی بهویژه سرما و بجیتنان بسیاز سای انن؛ بابراین جنانی د 
دمای محیط در مرحله شکوفه‌دهی کمتر از مقدار بحرانی گردد» بسته 
محصول دچار صدمه شده و گاهی خسارت ناشی از سرما اقتصاد 
کشاورزی یک منطقه را به‌شدت تحت تأثیر قرار می‌دهد (66 121 
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7 .01 
در ایران اگرچه آمار دقیقی از میزان خسارت سرمازدگی در 
سال‌های مختلف در دست نیست ولی با تحلیل غرامت‌هایی که 
هرساله توسط صندوق بیمه کشاورزی پرداخت می‌شود می‌توان 
برآوردی از حداقل زیان‌های وارده به کشاورزی کشور از این رهگذر را 

به دست آورد. از بین ۱۲ محصول استراترژ 
بوده‌اند بیشترین خسارت مربوط به سرمازدگی است که به‌طور متوسط 
۸ درصد کل غرامت پرداختی بوده است (2014 ,1621111). بر 
اساس گزارش معاونت باغبانی در اثر پدیده سرمازدگی در سال ۹۶ 
بالغ‌بر ۷۱۰۰ میلیارد تومان (در سطح ۶۶۱ هزار هکتار) و در سال ٩۷‏ 
بیش از ۱۱۰۰۰ میلیارد تومان خسارت (در سطح ۶۳۵ هزار هکتار) به 
باغات کشور وارد گردیده است به‌طوری‌که ۴۷۳۰ تن محصول باغی 
کشور از بین رفته است ( طمتاممتممطع۱۷ هنوخ 
8 اجهن) اجعجطع۲(۵۷۵1۵). در سال ۱۳۹۸ در ایران بیشترین 
میزان تولید در ین محصولات باغی مربوط به سیب با تولید حدود ۴ 
میلیون تن و سهم ۱۷/۱ درصد از کل میزان تولید محصولات باغی 
بوده که استان آذربایجان غربی با سهم ۲۷/۵ درصدی در رتبه اول 
تولیدکنندگان سیب کشور قرار داشته است (2020 .۵1 6 نحتطش). 
میزان خسارت سرمازدگی در سال ۹۷ در استان آذربایجان غربی 
۱ میلیارد ریال (در سطح ۱۱ هزار و ۳۹۰ هکتار) بوده است 
(2018 مهن اصمم‌ماهع۱ صمتاهمتصمطهع۱۷ اسان تتهض) 
بنابراین پرداختن به موضوع سرمازدگی محصولات کشاورزی و ارائه 
راهکاری برای کاهش خسارت‌های ناشی از آن اهمیت زیادی دارد. 

به دلیل اهمیت موضوع در پژوهش‌هایی که به‌عنوان مثال توسط 
امیدوار و دهقان بنادکی ( ,3202001 صقطع‌دهن1 فصه هنن 
2 و حصاری و همکاران (2015 ,.01 61 271:ع۳) صورت گرفته 
به تحلیل ویژگی‌ها و انواع یخبندان در ارتباط با محصولات کشاورزی 
پرداخته‌شده و یخبندان‌ها از لحاظ شدت در ارتباط با محصولات 
کشاورزی دسته‌بندی و بر مبنای آن تحلیل مکانی - زمانی سرمازدگی 
انجام‌شده و نقشه‌های وقوع سرمازدگی برای مناطق موردمطالعه 
تهیه‌شده است. با توجه به آهمیت پیش‌بینی تاریخ‌های فنولوژیکی 
برای پیشگیری از خسارات ناشی از متغیرهای هواء مدل‌سازی 
پیش‌بینی زمان وقوع مراحل فنولوژی محصول از اقداماتی است که 
موردتوجه قرارگرفته است و از روش‌های مختلفی مانند روش شبکه 
عصبی مصنوعی (2010 ,.۵1 ۶۶ طقصدم‌صد ۷۵2 ( ,.۵ ۵1 5212 
4 و روش تشخیص مراحل فنولوژی محصولات از طریق تصاوبر 
ماهواره‌ای و شاخص‌های پوشش گیاهی (2013 ,./۵ 6 تتاطل)» 
(2018 .اه 6۶ ۸۱۷۵۲62 62ج ۷۵۸) استفاده شده است. 
پیش‌بینی تاریخ دقیق رخداد گلدهی نیز نقش مهمی در کاهش 
خسارت ناشی از دماهای پایین دیررس بهاره دارد. تاکنون محققان 
مدل‌هایی مانند 1 32 (2011 ,۵ 6 فلدسه‌امنهع) و 


دینامیکی (2017 ,تلحفصطه مه تحطلفدن طعاله۳) را برای 
پیش‌بینی زمان گلدهی ارائه داده‌اند. هور و آهن ( معط 4ص ۲۲ 
15 نی با تفیق دو مدل فنولوژی و هواشناسی (مدل "0073 
تاریخ اولین گلاهی گیلاس» هلو و گلابی را در کره جنوبی پیش‌بینی 
کردند. نتایج تحقیقات آن‌ها نشان داد که روش مذکور نمی‌تواند تاریخ 
شروع گلدهی را در منطقه موردمطالعه به‌درستی پیش‌بینی کند. 

در مطالعات دقیق‌ت علاوه بر تعیین بهترین مدل برای پیش‌بینی 
گلاهی, تاریخ احتمال وقوع آخرین یخبندان نیز برآورد شده است 
(2007 ,له ۵ تلقصفکل» (2010 وله ۶۶ طعه2هتد۳) 
(2016 ۹۵022/۷72 0ص 2129۳601 ۷). بررسی روند تغییرات 
درجه حرارت در بلندمدت» روشی است که برای پی بردن به تغییرات 
ایجادشده در تاریخ وقوع گلدهی و نیز تغییرات در تاریخ وقوع آخرین 
یخبندان توسط مارتینز - لوسچر (2017 ,567 نا [-67ط1] ۷۲2۲ و 
ویتاسه و همکاران (2018 ,.۵1 2 ۷1]2956) استفاده‌شده است. 
نصر اصفهانی و یزدان پناه ( مطقصدم‌صه۷۵20 0صد تصحطح/ع۴ مدا 
9 دریافتند که می‌توان با کمک مدل فنولوژی و شبیه‌سازی دما 
به‌وسیله مدل ۰۷۷۲۳ هشدار سریع ۴۸ ساعته رخداد سرمازدگی در 
باغات منطقه موردمطالعه را با دقت کافی انجام داد. هدف اصلی و 
نهایی همه تحقیقاتی که در زمینه سرمازدگی صورت گرفته کم کردن 
خسارات ناشی از سرمازدگی است این تحقیق نیز در راستای همین 
هدف انجام شده است. باوجود پژوهش‌های زیادی که درزمینه 
سرمازدگی به محصولات به‌ویژه در زمان گلاهی صورت گرفته است 
اما همچنان لزوم ایجاد یک سیستم هشدار سریع وجود دارد تا برای 
هر باغ بتوان پیش‌آگاهی را انجام داد پس چون می‌خواهیم 
پیش‌آگاهی از سه روز قبل انجام شود یکی از اهداف این پژوهش 
ارزیابی دقت و صحت شبیه‌سازی‌های ۷۲ ساعته دما به‌وسیله مدل 
8 است. از طرف دیگر بایدبدانیم در هر یاغ محصول در کدام 
مرحله فنولوژیکی قرار دارد تا بر اساس دمای بحرانی هر مرحله اعلان 
خطر سرمازدگی به کشاورز انجام شود. بنابراین دومین هدف این 
پژوهش تشخیص مراحل فنولوژی از طریق تصاویر ماهواره‌ای لندست 
۷ ۸است. 


مواد و روش‌ها 

مناطق مطالعاتی 

منطقه مطالعاتی در ارزیابی مدل ۲۷1۳۲ 

(کوچک‌ترین شبکه مدل) مورد بررسی فرار می‌دهد. این منطقه شامل 
استان‌های آذربایجان شرقی» آذربایجان غربی» اردبیل و زنجان است. 
در این مطالعه از ۵۱ ایستگاه هواشناسی سینوپتیک واقع در شمال 
غرب ایران استفاده تلم 


ابن حجازی و همکاران. پیش آگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۲۷۷1۸۲ ۱۳۱ 


الف) موقعیت جغرافیایی کوچک‌ترین شبکه مدل 
ام مه عمط ن ممتایه۱ ای‌نطام‌منهمع6 (2 
0رع 5 ۱۵06 


ب) شبکه‌های محاسباتی انتخابی برای پیش‌بینی عددی وضع هوا 
مه ماه تما عل‌زرم امصمتاه) ری امعم (رط 
(۱۲۷۲۳) ۳۲۳60:6600 


ج) موقعیت ایستگاه‌های سینوپتیک 
»م۱ عصمتاهاه عتاممصک رم 


2190۳0۳۷ 


20۳0۳0۷ 


۳۳ ۱۱ 
60163 6 6, ۶ ۱ 


1 
۹-۹99۶ 4۹99090 تفشاف اتف فا 


0 س جه 10 
د) مدل رقومی ارتفاع 
(۲۳/۱۷) [ع۱۷۲۵0 م۵ز)۲۱۵۷۵ آه)نعز۱ ر0 


۱۷ 


۳۳ 


انش شرفت 


۲ اف تا 
شکل ۱- الف) موقعیت جغرافیایی کوچک‌ترین شبکه محاسباتی - شمال غرب ایران- که ارزیابی شبیه‌سازی مدل "۱12 در محدوده ان شبکه 
انجام شد. ب) شبکه‌های محاسباتی مدل ۰۷۷1۲ ج) موقعیت ایستگاه‌های سینویتیک که برای ارزیابی شبیه‌سازی‌های مدل ۲۷1۸۲ استفاده شدند. 


د) مدل رقومی ار تفاع شمال غرب ایران (0۱05://25.1) 
)۱۵۲۳۱۷6۵ متام احصمناه) نمی )عماامصصی معط ۵۶ جمتایه۱ انامه و۵عع (۵ وم ۱۵۲ ۲۷۷ ۵ م0ام امصماه)نم‌صم (ه .1 مساهز۲] 
۶ ۱,۵۵0 (0 راهم میاه امانونل مه اه )مها رن بان ماه ۱۵۵0۵۵ ۲۷۷۱۲۸۲ ۵ اممجرووم‌ووه )5ه06صنط طمزظ۳ ۸ مضه 
۲۵۹۰//۵16۰1۴(۰) (0۵معه ۲۷۷۱۸۲ ۵ نی دعمععه ما معا ۱۷۲۵۲۵ عممتاهای »ا۳۵۵۵ه 


شکل‌های ۱ الف» ۱ ج و ۱ د موقعیت جغرافیایی» موقعیت شبکه 
ایستگاه‌های سینوپتیک و نقشه مدل رقومی ارتفاع شمال غرب ایران 
مورداستفاده برای ارزیابی مدل را نشان می‌دهد. 


مقادیر شاخص پوشش گیاهی برای یک پیکسل ۳۰ متری واقع 


در مرکز باغ آزمایشی سیب در اداره تحقیقات کشاورزی کهریز استان 
آفربایان خزی فر فمال غرب‌ارآن 3930:3000 3۳ 
(57272 برای سال‌های ۲۰۰۷-۲۰۱۶ استخراج شد (سکل ۲). بر 
اساس آمار بلندمدت ایستگاه تحقیقات هواشناسی کشاورزی کهریز 
ارومیه طی سال‌های ۱۳۶۵-۱۳۹۵ میانگین بيشینه دما ۱۷/۴ درجه 
سلسیوس میانگین حداقل دما ۷/۲ درجه سلسیوس, میانگین رطوبت 
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نسبی ۵۵ درصد و مجموع بارش سالانه ۲۸۹/۶ میلی‌متر است 
(صاط صمنماه لصمناماه لصا اع ۰227 ۰)۳۱)0://۷۷۷۷ نوع خاک آن 
سیلتی لومی است. این منطقه یک ناحبه کوهستانی با فرورفتگی 
وسیع و شیب ملایم است. تراکم باغ سیب رقم مالینگ ۸۲۳ اصله در 
هکتار (فاصله ردیف‌های کاشت ۴ متر و فاصله درختان روی ردیف ۳ 
متر) و تاریخ کاشت درختان در سال ۱۳۸۰ بوده است (مرکز تحقیقات 
هواشناسی کشاورزی کهریز ارومیه» ۱۳۹۵). 


داده‌ها 
داده‌ها و مدل ۲۷۷1۵۲ 

مدل پیش‌بینی و تحقیق وضع هوا (0۷1517) نسل بعدی سیستم 
پیش‌بینی عددی میان‌مقیاس آب و هوایی برای کاربردهای تحقیق و 
پیش‌بینی جوی است. ۱۷15۲ می‌تواند شبیه‌سازی‌ها را بر مبنای 
شرایط جوی واقعی (به‌عبارت‌دیگر بر پایه مشاهدات و تحلیل‌ها) و 
شرایطی که به‌صورت ایده آل و دلخواه ایجادشده است. تولید کند 
ینامیا 

ایق فطالعه قوانایی قبیاسازی دمای حداقل به‌فسله متل ۲ 
را در شمال غرب ایران مورد بررسی قرار داد. با توجه به هدف 
پژوهش که پیش‌آگاهی سرمازدگی است. ارزیابی پیش‌یابی مدل برای 
دمای ۲ متری در ساعت )۲71 ۳ انجام گرفت. به دلیل اينکه درختان 
سیب در فصل بهار و در طول دوره گلدهی در معرض دماهای پایین و 
خطر یخبندان قرار دارند شبیه‌سازی ۲ روزه (۷۲ ساعته) دمابامدل 
۷۷۳۷ نسخه ۴ برای ۱٩‏ روز نمونه از یک هفته قبل از 
شروع دوره گلدهی از ساعت ۱۲ روز ۲۸ مارس تا ۱۶ آوریل ۲۰۱۲ 
انجام شد. ۱۲ ساعت اول برای زمان تنظیم (015-0۳ه) در نظر گرفته 
شد. شبکه محاسباتی برای شبیه‌سازی دمای ۲ متری به‌وسیله مدل 


۰ ۳۳۳۳۲ بر سس 


«(5 


میس ۰ ۲ 


8 تب ۳ 
1-۳ سس 


عددی ۰۷۷۳۳ متشکل از سه شبکه تودرتو (لانه‌ای) است (شکل ۱ 
الف) که مشخصات آن به شرح زیر است: 

شبکه بزرگ دارای تفکیک افقی ۲۷ کیلومتر است و محدوده 
۲۰-۰ درجه شمالی و ۲۰-۶۱ درجه شرقی را پوشش می‌دهد. تعداد 
نقاط این شبکه ۱۶۷۱۲۲ نقطه است که مساحتی بالغ‌بر ۳۲۹۴۲ 
۲ را شامل می‌شود. شبکه متوسط دارای تفکیک افقی ٩‏ 
کیلومتر است و محدوده ۲۶-۴۲ درجه شمالی و ۲۶-۵۵ درجه شرقی 
را پوشش می‌دهد. تعداد نقاط این شبکه ۲۳۲۰۱۹۶ نقطه است که 
مساحتی بالغ‌بر 160 ۱۷۶۴ ۲۰۸۸۲ را فرامی‌گیرد. شبکه کوچک 
دارای تفکیک افقی ۲ کیلومتر است و محدوده ۲۵-۴۰ درجه شمالی و 
۳۳-۰ درجه شرقی را پوشش می‌دهد. تعداد نقاط این شبکه 
۸۱ نقطه است که مساحت آن بالغ‌بر ۲۵۴۹ ۶۶۰1 است. 
تفکیک افقی داده‌های ناهمواری و کاربری زمین برابر ۳۰ ثانیه (حدود 
۱ کیلومتر) است. شرایط اولیه و مرزی از داده‌های سامانه مدل‌سازی 
تمام کره‌ای موسوم به 015 در بازه‌های ۳ ساعته و با تفکیک افقی 
نیم درجه از مرکز ملی داده‌های محیطی (0۲1) گرفته شد. مدل 
5 (ععاور5 ]25»ع۳۵۳ 010021 ۲6)» یک مدل پیش‌بینی است 
که توسط مراکز ملی پیش‌بینی محیطی (02۳) تهیه می‌شود. تعداد 
زیادی از متفیرهای جوی و زمین - خاک در این مجموعه داده موجود 
است. از دماهاء بادها و بارش تا خاک رطوبت و غلظت ازن اتمسفری 
(10۸۸). بر پایه پژوهش‌های پیشین پیش‌بینی عددی وضع هوا در 
ایران (2013 ,./ه 2 ت0همز20۱) ( ,تلقه2ظ 40صده صحتعگعطت 
3 (2014 ,اه 6 تع 22 (2012 .اه ۶ 220ه) از 
طرح‌واره‌های فیزیکی پیکربندی 16۷۲۷1 برای مدل ۱۷۲۳ در این 
تحقیق استفاده شد (جدول ۱). 


شکل ۲- تصویر سمت چپ: باغ آزمایشی سیب اداره تحقیقات هواشناسی کنساورزی کهریز استان آذربایجان غربی منبع تصویر: 12۲10 200210 
تصویر سمت راست موقعیت جغرافیایی باغ مطالعاتی در نقشه ایران 


۵0۱6 ۷۷ و۸ ص۳۵۷ صوزنه۸2۵0 ۲۷۲۷۲۵۵۲ ر2تتطاهک ۵ فصن آماممهمن مامروره :۱۵ عظ) م0 -2 عم۲اع۲1 


۵۵ ۲۱۳۵۱ 0تفصمرن اه)صمصتمه موی مامرمره معط ۵۶ جمتای‌م۱ امتصام‌موممع )صعز۱ مط) جر نو مزهممع بععتمو قفا بماممن 


اين حجازی و همکاران. پی شآگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۷۷1۸۳ ۱۳۳ 


جدول ۱- طرح‌واره‌های فیزیکی انتخاب‌شده جهت مدل‌سازی با مدل ۲۷۷۲۸۲ 
۲ ۳ هام۵0 م۲ عمصمطامی هدیدان 5۵16060 -1 12016 
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با استفاده از پیکربندی منتخب (جدول ۱) شبیه‌سازی دمای ۲ 
متری بامدل ۷۷۲ برای ۷۲ ساعت (۳ روز) انجام و مقادیر 
به‌دست‌آمده در کوچک‌ترین شبکه در ساعت ۲7۲0 ۳ برای هر یک از 
ایسشگاههای. هیاشتانی بستوشک قمال غربه اپران.- کته داوه‌های 
آن از سازمان هواشناسی ایران (11۷10) اخذ شد-در نرم‌افزار متلب 
با تابع 106703 و با روش 110627 درون‌یابی شد (معادله۱). البته 
روش‌های دیگری نیز برای درون‌یابی مورد ارزیابی قرار گرفت که 


روش مذکور بهترین عملکرد را داشت. 
)۱ (0مصامصر0 0 ۷ ,20 ۷ ور ۷ ,)1010103 عون ۱۷ 


که ۲ عرض جفرافیایی» ۷ طول جنرافیایی؛ 2 زمان و ۷ 
متغیری است که درون‌یابی برای آن انجام می‌شود و در این مطالعه 
منظور دمای ۲ متری است. 264 و ۷۹ طول و عرض جغرافیایی مکان 
و 24 زمان موردنظر برای درون‌یابی متغیر است. به‌منظور ارزیابی 
عملکرد مدل از مفاهیم آماری مانند میانگین قدر مطلق خطا 
(۱1۸2), خطای میانگین مربعات (1155» مجذور خطای میانگین 
مربعات (18۷]5۳)» شاخص امتیاز مهارتی (0۷1555» شاخص میانگین 
اریبی (۷18» ضریب همبستگی (*) استفاده‌شده که با استفاده از 


فیزیک مدل 
مدز ایزوود۳ 
طرح‌واره میکرو فیزیک 
۷۱۱۵0۵0( 
طرح‌واره تابش طول موج بلند 
0۵[ 
طرح‌واره تابش طول موج کوتاه 
۱0۳۵۹۵۹0 
طرح‌واره لایه سطحی 
۱ 
طرح‌واره سطح زمین 
۱9۰۱۱۱۵۱۱ 
طرح‌واره لایه مرزی سیاره‌ای 
۲ ۵021-12 )م۳۱ 
طرح‌واره همرفت 
مرن 
گام زمانی 
سم 11۳06 
تعداد ترازهای قائم 
عاه۱6۷ ۷۵۲۵۵۵ ]0 ۱0 
طبقه‌بندی کاربری اراضی 
60۲۵۲ 2110[ 


معادله‌های ۲ تا ۷ تعریف می‌شوند: 
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در معادله‌های فوق :) مقادیر دیده‌بانی شده و :۲ مقادیر پیش‌یابی 
شده و 7 تعداد داده‌ها است.0 میانگین مقادیر مشاهداتی و۲ 
میانگین مقادیر شبیه‌سازی‌شده است. 

۷۸ ۷۲5۳۲ ۷۲۹۳۴ و 1۷555 هرکدام به شیوه‌ای مقدار 
اختلاف بین شبیه‌سازی مدل و مشاهدات را نشان می‌دهند و هرچه 
مقادیر کمتر باشد یعنی کارایی مدل بیشتر است. شاخص 1۷8۳ 


۴ نشریه آب و خاک. جلد ۰۳۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


علاوه بر دارا بیدن علامت مثبت (بیش برآوردی مدل) یا منفی (کم 
برآوردی مدل) مقدار انحراف از مقادیر دیده‌بانی شده را نیز ان 
می‌دهد (2019 ,.1 6 تصدطد1ه28-۲). *1 معادله خطی بین 
شبیه‌سازی مدل و داده مشاهداتی را نشان می‌دهد ( ۶ تصدامطان 
7 ,./). مقادیر ضریب همبستگی همواره عددی بین ۱ و ۱- 
است. هرقدر ضریب همبستگی به ۱ نزدیک‌تر باشد عملکرد مدل در 
شبیه‌سازی تغییرات پارامتر بهتر است. 


داده‌ها و ارزیابی تصاویر لندست 

در اين مطالعه از آمار فنولوژی - زمان شروع و پایان مراحل رشد 
- که به‌صورت سالانه در باغ سیب اداره تحقیقات هواشناسی 
کشاورزی به ثبت رسیده (جدول ۲) استفاده شده است. 

به‌علاوه ۷۷ تصویر ماهواره‌ای لندست ۷ سنجنده +۳۲1۷ و ۴۱ 
تصویر لندست م۸ سنجنده لام از سایت 
۲ ۵ در ماه‌های مارس تا سپتامبر 
سال‌های ۲۰۰۷-۲۰۱۶ (مطابق با مراحل فنولوژی سیب در ایستگاه 
کهریز) با مشخصات شماره گذر و ردیف به ترتیب ۱۶۹ و ۲۴ و قدرت 
تفکیک مکانی ۲۰ متر دانلود شد. در این تحقیق پیش‌پردازش‌هایی 
شامل تصحیح اتمسفری و رادیو متریک با روش ۳1۸۸۹۲1 بر روی 

محاسبه شاخص پوشش گیاهی از طریق شاخص ۱۲2۷۲ انجام 
سطحی با توجه به طول موج تعمیم داده شوند. این اشتقاق تابع 
خصوصیات نوری برگ‌ها و ذرات خاک است ( ,۵1 6 ۱۷۲۷۲6 
5 علاثئم بازتاب طیفی. اطلاعاتی را در مورد وضعیت ترکیبات 
بیوژئوشیمیایی و ساختار برگ و تاج پوشش گیاهی آشکار می‌سازد 
(2012 ,عاه]۲). 


شاخص نرمال شده تغییرات پوشش گیاهی (۲«۷) (معادله )٩‏ 
پرکاربردترین شاخص در مطالعات مربوط به پوشش گیاهی است 
۱۱ 
8 (2013 ,صمتاهاطاعدن لمع 122 ماعطحق) ( رعتعاعع۱۷]0۳ 
6 شاخض 3۳0۷/۲ بر اساس سبت باندی بین باندهای ماذون 
قرمز و قرمز محاسبه می‌گردد و دامنه تغییرات آن بین یک و منفی 
یک است (2016 .۵ 6 ۲۵۳10۲01 [۳۵2۵). فرمول ذیل برای 
محاسبه ۲۷۲ استفاده می‌شود (۲۰). 

)0 (ومی؟ + مر /(ومی؟ - ور < ۸۷۷ 

* نشان‌دهن ده انعکاس در هر ناحیه مشخص از طبیف‌های 
جمع‌آوری‌شده است (2013 ,.۵1 6 626۲). 

در این تحقیق شاخص ]۱۲2۷ با استفاده از داده‌های لندست ۷ 
سنجنده +1۷1 و لندست ۸ سنجنده 0 در طی دوره رشد درخت 
سیب در ماه‌های مارس تا سیتامبر سال‌های ۲۰۰۷-۲۰۱۶ برای 
منطقه مطالعاتی در نرم‌افزار 1۲۷ محاسبه شد. 


ارزیابی پیش‌یابی دما 

با بررسی شاخص‌های ارزیابی مدل ۷۷۴۲ چنین نتیجه‌گیری 
می‌شود که میزان صحت شبیه‌سازی‌های ۰۲۳ ۳۸ و ۷۲ ساعته 
درمجموغ ایستگاه‌های منطقه شمال غرب ایران تفاوت زیادی یاهم 
ندارند ولی به‌طور کلی ازنظر صحت. شبیه‌سازی ۲۴ ساعته در رتبه 
اول و شبیه‌سازی‌های ۴۸ و ۷۲ ساعته در رتبه‌های بعدی قرار دارند 
(جدول ۳). شکل ۳ نشان می‌دهد هر چه زمان شبیه‌سازی افزایش 
موزنید پراکندگی تقاط حول خط رگزسیون پیشتر من شود 


جدول ۲- زمان وقوع مراحل فنولوزیکی باغ سیب ایستگاه تحقیقات هواشناسی کشاورزی کهریز بر اساس داده‌های زمینی 
2 ۱6۱ و 60و صمتاهاه همع مزع ما0۲0ع)مص0 ۸۵ 1۵۱2 ۵۴ فصن مامرهرد ۵۶ عععهاه ام‌زعهامصمدام 01 11۳6 -2 12016 


روز سال 
(1۵۲) ۲۵۵۱۲ 0۴ 1۱۵۸6 سال 
رسیدن میوه رشد میوه گل دادن شکفتن جوانه تشکیل جوانه وه 
مصنصه مرن ات مهم ات مصته۲ ۲۱۱۵ هصنصهمععماظ .. صمت)ههم 0اظ 
251 135 12 109 95 20007 
236 10۳4 91 90 83 20058 
257 124 103 96 59 200 
258 118 91 / 609 2010 
259 125 115 103 94 2011 
251 132 115 106 91 2012 
259 10۳4 96 84 69 2013 
254 114 91 90 69 2014 
273 118 100 92 13 2015 
272 121 107 93 82 2016 
257 120 105 94 82 میانگین (۱۷۲6۵0) 


ابن حجازی و همکاران. پیش آگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۲۷۷1۵۳۲ 


۱۳۵ 


جدول ۳- نتایج شاخص‌های مورداستفاده در ارزیابی شبیه‌سازی مدل ۲۷۲1۸۲ 
عجرم میتی امه 1۲ ۲۷ ۵۴ عصتا ماو صا 0مکده ومعلصا ۵ دااناععع1 -3 12016 


۳( و۷ ۳۱۷۲9۲ ۷۲9۲ 

1-3 ۵4 0۳ 0909 

11-385 2392 ۱ 54 11617 

1211 005568 ۰-35" 1-۱ 
کمماهل ااد 


085 < 3.122 + 05216 ( 


085 - 2.811 + 0.5465 ۵ 


مدت شبیه‌سازی (00-۵6100 صمتادلبصنق) 
248 ۴ ساعته (24-1000) 

282 ۸ ساعته (48-000) 

288 ۳ ساعته (72-00) 
جهتاهء ااد 5وزاهاو اج 


085 2 2.861 + 0.5457 ۵4 


شکل ۲- معادله و نمودار خط ر گرسیون داده‌های مدل ۷۷3۲ و داده‌های مشاهده‌شده در منطقه شمال غرب ایران در شبیه‌سازی‌های ۰۲۴ ۴۸ و ۷۳ 


ساعته 


۱ 
کصمز)هاصصتی عمط 72 200 48 م24 از ما۲ مه رز مصمتاهای »موی ص۲۵0 عوهای ماه ماهنت مرها 


جدول ۴- خلاصه مدل 
۱۷۲۵061 -4 12016 


۵ 0۱۵6 0۲ ۳۸۲۵۲ :۹۲۵ 5002۲6 1 0ع)وازل ۸ 
خطای استاندارد برآورد ضریب تعیین تعدیل‌شده 
2-2-1255 0345 
3206821 0359 
11100--2- 0305 


۷۲۵0۵ ۵۵۵۲ ۳-۷ 2۱06 ۰ 6 

ضریب تعیین مقدار-۳ خلاصه مدل 
316 0۳ 0 با _-323932 3«( ۷4 
0359 0 . اووو.0 ۷46 
306« 0.- 0.553 ۷2 


: (پیش‌بینی کننده‌ها) ۳۲۵01000۲5 .6 48 ۷۴ ,(ثابت) (08۵/2۳0)) :(پیش‌بینی کننده‌ها) ۳۲601060۲5 .0 24 ۷۲۶ ,(ثابت) (8]2ظ0)) :(پیش‌بینی کننده‌ها) ۲۲6۵0160015 .2 
2 ۷ , (ثابت) اصداعدمن)) 


میزان ضریب همبستگی و ضریب تعیین بین متغیرهای مقادیر 
دمای حداقل به‌دست‌آمده از مدل ۷۷1۵۳ در کل ایستگاه‌ها و مقادیر 
مشاهده‌شده در جدول ۴ نشان داده‌شده است. منطقه شمال غرب 
ایران یک منطقه کوهستانی است و مدل ۷۷۳ در مناطق کوهستانی 
در شبیه‌سازی کارایی کمتری دارد این امر باعث شده تا درمحموع 
مقدار ضریب همبستگی در روزهای موردمطالعه در حد بالایی نباشد. 
با این حال معیار معنی‌دار بودن يا نبودن میزان همبستگی, مقدار 
احتمال آن است که باید با سطح خطای ۰/۰۵ مقایسه شود با توجه به 
اینکه در جدول ۴ مقدار ۳-۷۵10 از ۰/۰۵ کمتر است درنتیحه میزان 
همبستگی شبیه‌سازی‌های ۲۴ ۴۸ و ۷۲ ساعته با دمای مشاهداتی 
معنی‌دار است. همچنین در جدول ۵ با توجه به اينکه مقدار ۳-۷106 
یا همان 9121110206 کمتر از ۰/۰۵ است می‌توان گفت مدل 


رگرسیونی ارائه شده معنی‌دار است و مدل به‌خوبی قادر به شبیه‌سازی 


تغییرات متغیر وابسته است. 

جدول ۶ ضرورت وجود هر یک از متفیرها و به عبارت دیگر 
معنی‌داری ضرایب متفیرهای مستقل در مدل‌های رگرسیونی برای 
شبیه‌سازی‌های مختلف را نشان می‌دهد. همان گونه که در ستون 
1 جدول مذکور مشاهده می‌شود ضریب ثابت 0صداععهع) 
و ضریب متغیر هر دو برای پیش‌بینی متغیرهای وابسته معنادار هستند 
و می‌توان از آن‌ها در معادله رگرسیون استفاده کرد. 

با توجه به اینکه هدف اصلی در این تحقیق ارزیابی شبیه‌سازی 
دمای حداقل ۷۲ ساعته است. توزیع جغرافیایی شاخص‌های 
محاسبه‌شده برای ارزیابی شبیه‌سازی متغیر مذکور در ایستگاه‌های 
منطقه شمال غرب ایران در شکل ۴ نشان داده‌شده است. الگوی 
نقشه شاخص‌های خطا مطابق با کوهستانی بودن منطقه است و با 
تحلیل نقشه‌ها به‌صورت استانی این نتایج به دست آمد: 


۶ نشریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


جدول ۵- تحلیل واریانس 
۸۸ ۸0۲۷ -5 12016 


12۰ 7 6 ۷۲6۵7 
سطح معنی‌داری میانگین مربعات 

+302 520639 00000 
9115 

229.033 2-1 00000 
«14 

"۹-111 (29 0000۹ 
968 


500۵ 0 ند 
۳ ع مربعات ۸ ۷۲0061 
۱ ۳ مدل تحلیل واریان 
پراکندگی وارباس 
4 ۲۷۷۲۲ 168۳625107 
1 ۹4502 4 ۱۷12 رگرسیون 
4 ۷۷۲۲ 165100721 
959 9313503 4 ۷78 باقیمانده 
4 ۷۷۸۲۳ 1021 
960 119.05 4 1771817 کل 
8 ۷۷۸۲ 162655100 
1 29033 48 1۷1۵ رگرسیون 
8 ۷۲۷۱۲ 165101021 
960 903373 48 ۱۷131 باقیمانده 
۷۷۸۲ 1021 
961 ۱۳۹۹ 48 1۷۵ کل 
2 ۷۷۲ 168۲62510۳7 
1 9.111( 2 ۱۷1 رگرسیون 
2 ۷۲ 1651010721 
960 9291-183 2 17718 باقیمانده 
2 ۷۷۸۳۲ 1021 
961 1-94 2 ۷7 کل 


:(پیش‌بینی کننده‌ها) ۳۲۵010/0۲5 6۰ 24 ۷1۴ ,(ثابت) (00۵18۳0)) :(پیش‌بینی کننده‌ها) ۳۳۵016۲0۲5 .0 (مشاهده‌شده) ۵096۲۷6۵ :(متغیر مستقل) ۷۵11۸016 1۵۵۵0601 .2 
2 ۲ , ( ثابت) (احماعطمت) :۳:606)0۲9 :0 48 ۷۲۶ , (ثابت) (اصفاعجمت) 


جدول ۶-جدول ضرایب * 
«ماممون: 00011 -6 12016 


عاصعتن 00115 02۲01760ظه5 
(ضرایب استانداردشده) 


.12 
سطح 1 ۳2 
معنی‌داری 
0.000 12619 
0.000 2.13 0388 
0.000 123270 
0.000 23.205 0599 
0000 12461 
0.000 20.75 03553 


عاصهنن]1ه0) 2۵027-01760 عصر] 


(ضرایب استاندارد ننسده) 


ی عاصعهت 011 1۷۲۵061 
۳ نت 

خطای استاندارد ۳ 
02217 2-1 (ثابت) (اصعاعصم۲])) 
0024 036 ۳/۳۳۹4 
00229 2-11 (احعاعصیت) 
94«( 347( ۱۵ 
002351 2-122 (احعاعصیت) 
00025 0322 ۷۷۰۲۰2 


(مشاهده‌شده) 00۹6۲۷6۵ :(متغیر مستقل) ۷۵112016 1۵۵۵0601 .2 


ناهمواری‌های عمده در استان آذربایجان غربی در قسمت غربی 
آن یعنی در نواحی مرزی ایران با ترکیه و عراق واقع گردیده است. 
رتفاعات عمده این ناحیه را کوه‌های مرزی ماکوه قطور, خوی در 
شمال استان و کوه‌های غرب ارومیه. اشنویه. پیرانشهر و سردشت در 
جنوب استان تشکیل می‌دهند ( ,طمتا2012ع0۳ )۳06 0ص ص۳۱ 
8 به سمت نواحی کوهستانی مذکور همبستگی اغلب به کمتر از 
۷ و خطای ۷۲۸۲ مدل به بیش از ۲/۵ درجه سلسیوس میرسد. در 
استان آذربایجان شرقی به سمت کوه‌های قره داغ در شمال تا سبلان 


در شرق و سهند در مرکز از میزان همبستگی کاسته شده و خطای 
شامازی تالا هی رود شکور آنیشگاه لیا کد سم رنه 
رشته کوهستانی قره داغ در آن قرار دارد میزان همبستگی به ۰/۶۴ و 
میزان خطای ۷۲۸۳ مدل به بیش از ۳/۵ درجه سلسیوس میرسد. در 
تام ارف کر کنر آنسشگاه‌ها ریب ,هسسگی تال هو 
خطای ۷1۸۳ مدل تا ۲/۵ درجه سلسیوس است. در استان زنجان در 
ایستگاه‌های بین کوه‌های طارم و سلطانیه خطای ۷1۸ مدل بالای 
۵ درجه سلسیوس است. 


ابن حجازی و همکاران پیش آگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۲ ۵ ۵ ۱۳۷ 


ِِكِِ« ٍِِ۹ِ«آىكآكآ‌"«- 
الف) میانگین قدر مطلق خطا (۷۲۵) 
0۷۲۸۲ ۲۱۱۲۰۵۲ ع6امعوا۸ صهع]۱۷ (2 (10) «0تامامنه (ط 
«آ«/«. 4 تا تاه خی 


پ) خطای میانگین مربعات (01519) ت) مجذور خطای میانگین مربعات (111519) 
(۷۱۵۲) 0۲ هفوک صفه۷۲ ه (۳]61) 0۲ محمیو 5 صهع]۷ 13001 (0 


بت ند ای اف ۰ 


۰ 9 
ث) شاخص امتیاز مهارتی (۷۲۹55) ج) شاخص میانگین اریبی (117) 
ردوق۵ م66 الاک 0عسهنوک مهه ۵ (۷۱۳۵۲) ۲۲۵۲ ععظ «هع۱۷۲ ( 
شکل ۴- نقشه توزیع جغرافیایی شاخص‌های محاسبه‌شده در ارزیابی پیش‌یابی ۷۲ ساعته: الف) 2۷۲۸۰ ب) 35 پ) ۷۲57 ت) 381510 ث) 
۹ ج) ۷۲۶۳ 


(0 و(۷۲۸۳۲ 2 :صمت)م او )عق0صنظ سامط-72 ]۵06 ۱۵۲ ۲۷۷ معط مد 0ع)هانماه ومع10 مط) ۵۶ مرهصه ممتاطاتاعتل امومع -4 مساعز۲ 
۷5 ره و۱۷۲۳ (] ]1۷ (0 م۷۲ (0 مج 


۸ شریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


توزیع مقادیر خطای 1۷]9۳ در ایستگاه‌های مختلف» با خطای 
۸۲ هماهنگی دارد ولی به‌طورکلی مقدار خطای بیشتری را نسبت 
نشان می‌دهد. در ایستگاه‌هایی که خطای مدل کمتر باشد 
شاخص ۷595 به ۱ نزدیک‌تر است همچنین هر چه مقدار خطای 
مدل کمتر باشد شاخص 1198 به صفر نزدیک‌تر و درنتيجه پیش‌یایی 
واقعی‌تر است. 


ذشخیص مراحل فنولوژی باغ سیب از تصاویر ماهواره‌ای 
مقادیر شاخص ۱۷0۷1 برای یک پیکسل ۳۰ متری در مرکز باغ 
سیب ایستگاه هواشناسی کشاورزی ارومیه» مطابق با آمار مراحل 
فنولوژی ثبت‌شده در ایستگاه در سال‌های ۲۰۰۷-۲۰۱۶ و با استفاده 
از تصاویر ماهواره‌های لندست ۷ سنجنده +11۷ و لندست ۸ 


سنجنده آ1 به دست آمد. به‌طورکلی در منطقه موردمطالعه مقدار 
حین گل دادن به ۰/۴ و ۰/۵ و در دوره رشد میوه به بیشترین مقدار 
می‌رسد (جداول ۷ و ۸). با توجه به مقادیر ۱۷1 می‌توان نتیجه 
گرفت که هراحل فتولوژی ثبت‌شته از باغ سیپ ایستگاه کهریز با 
داده‌های شاخص پوشش گیاهی تطابق دارد. همان گونه که انتظار 
می‌رفت بعد از پیان دوره گلدهی با توجه به رشد برگ و نیز الا رفتن 
مقدار کلروفیل در ماه‌های گرم سال مقادیر بالاتری از شاخص 
11 را نسبت به دوره‌های خواب» شکفتن جوانه و گلدهی شاهد 


جدول ۷- مراحل فنولوژی سیب ثبت‌شده در باغ سیب ایستگاه مرکز تحقیقات هواشناسی کشاورزی کهریز ارومیه و شاخص پوشش گیاهی 70۷1 
به‌دست‌آمده از تصاویر +51۷1 7 )120052 
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۱۵8۵65 +۲,۱۱۷۲ 7 )1.2۱062 
۱( مراحل رشد 1 ۱( مراحل رشد ۳ ۱( ۲ ره 
تا / تا 
9۲) 20۳ 4 ۲) ۱ ۲) و / 0 ِِ 
(+1:1 7 معهاه ٍِ (+1:1 7 معهاه ِ (ب11 7 هه ۲۵۷۲۲۵ ِ 
گل دادر خواب شکفتن جواز 
050 ات 25( 0292 12 /۰<۵ ۵ 0304۵ ِِ 7 2 
عصته/۳۲۱۵۱۳ ۱9۹۸ 4( 
۳ 3 داد ۴ واه 
03۹6 1 1-15 0353 ِِ 2 . 0381 ِِ 0 
0۵ ۲۲۷۱۸ 10/۵ 10/۵۵ 
034 ۰ 015 ۰ 0.657 0 02 ۰.۰ 0.391 ان 7 72-2 
۵ ۲۱۲۱۱۱۸ 0۵ ۲۱۲۱۱۱۸ 8 ۲۱۲۱۱۱۸ 
0497 015 ۰ ۰ 0.690 0 2 .0.467 7 
۱۹ ۱۵ ۱۹۹۵ 
05۱ ۰ 015 "۰ 0.709 02 ۰ 0.627 ات 2-7 72 
08 ۲۱۲۱۱۱۸ 0۵ ۲۱۲۱۱۱۸ ۵ ۲۱۲۱۱۱۸ 
0405 0 015 ۰ 0743 ۰ 2 ۰ 0724 ۰ ۳ 
۱۹ ۱۵۵ ۱۹۵ 
039 با 015 ۰ 0718 ۰ 2 ۰-۰-2 ۵-2 02 ۰ 10-3-07 
۵ ۲۱۲۱۱۱۸ 0۵ ۲۱۲۱۱۱ 8 ۲۱۲۱۱۱ 
0415 ۰ 015 ۰ 0654 انیت 12 ۰ 0.626 ۰« 11-7 
08 ۲۱۲۱۱۱ ۱ ۱ 
ْ 0 : ۵ خواب 
0432 0 14-216 0084 ین 892012 0237 5.32008 
08 ۲۲۷۱۸ ۱۱0۵۵۵ ۲۱۲۱۱۱۲ ۱9۸ 
0544 ی 2.2016 0359 6.2013 0352 208 
۱۱ 0۵ ۲۱۲۱۱۱۸ 0۵ ۲۱۲۱۱۱۸ 
0499 ۰ 19۹2016 086 9 762013 0448 ِ 92008 
۱۵ ۱ ۱۹۵ 
0۹86 ۰ 6 ۰ ۰ 0728 ات 013 ۰ 0.536 ۰ 268 
۵ ۲۱۲۱۱۱ 0۵ ۲۱۲۱۱۱۸ 8 ۲۱۲۱۱۱ 
040 ۰ 26 091 اس 013 ۰ 0.572 و 2-8 
۱ ۱۹۵۵ ۱۹۹۵ 
0461 ۰ 6 0669 ۰ 013 ۰ 0598 ِِِ 11-8 


ری ۱۱۱۱۱۱۰ 


۲۱۲] 20۵ 


م20 ]ن۳ز۲ 


ابن حجازی و همکاران. پیش آگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۷۷۲1۵۳ ۱۳٩‏ 
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۶۰ نشریه آب و خاک. جلد ۰۳۶ شماره ۱. فروردین - اردیبهشت ۱۴۰۱ 
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اتذال مراحل رشد تاریخ انذال مراحل رشد تاریخ 
ها ۲ ۳0۳۲) ۱9۹6 ۹22 6۳0۳۲۱۵ 1۱(9۹6 
و اد ۳ از 
0355 ات 19۹۹-2015 0۹0 ِِ 23 
۱ ۳1008۵ 
۳ 5 ی 
008 1 215 0546 ِِِ 12-2013 
۱۹۵ 1۳0۵۵ 
0212 ی 215 0369 0 2-4-2013 
8 ۲۱۲۱۱۱ ۱۱ 
ها شم 5 2 03721 0 30213 
۱ ۱۱۹۹۱ 
0513 7 15 0357 ها 1-213 
8 ۲۱۲۱۱۱ ۱۱ 
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2 10-9-2016 0۹61 0 2-1-14 
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در تحقیقات مرتبط با سرمازدگی محصولات» بر اساس آمار و 
یخبندان برآورد شده است اما روشی که در این مقاله برای پیشگیری 


ری ۱۱۱۱۱۵۰ 


از خسارت سرمازدگی در منطقه ارزیابی شد پیش‌آگاهی به کشاورز از 
طریق شبیه‌سازی دمای کمینه برای سه روز آینده با استفاده از مدل 
۳ قبل از زمان جوانه‌زنی بود که درنتیجه آن کشاورز بتواند با 


ابن حجازی و همکاران. پی شآگاهی سرمازدگی بهاره باغات سیب با استفاده از مدل ۷۷1۸۳ ۱۴۱ 


توجه به دمای بحرانی» فرصت کافی برای مدیریت وضعیت پیش رو 
داشته باشد. همچنین در مطالعات پیشین, بیشترین تمرکز بر روی 
تشخیص و پیش‌بینی مراحل فنولوژی به‌ویژه گلاهی بوده که نقش 
مهمی در پیش‌آگاهی سرمازدگی دارد این در حالی است که در ایران 
ایستگاه‌های هواشناسی کشاورزی اندک است و تاریخ دقیق مراحل 
فنولوژی محصولات مختلف در همه باغات ثبت نشده و درنتیحه 
آماری وجود ندارد تا بر مبنای آن مدلی برای پیش‌بینی گلدهی و سایر 
مراحل فنولوژی ارائه شود لذا در این تحقیق استفاده از تصاویر 
ماهواره‌ای برای تشسخیص مراحل فنولوژی سیب در ایستگاه 
هواشناسی کشاورزی کهریز ارومیه مورد بررسی قرار گرفت تا بتوان از 
تصاویر ماهواره‌ای برای تشخیص مراحل فنولوژی در مناطق فاقد 
ایستگاه هواشناسی کشاورزی در محدوده شهرستان ارومیه استفاده 
کرد. 

بنابراین این مقاله به بررسی امکان به کار بردن شبیه‌سازی ۲ 
روزه دمای حداقل با استفاده از مدل ۷۷۲۹۳ و نیز تصاویر ماهواره‌ای 
لندست برای پیش‌آگاهی سرمازدگی دیررس بهاره باغات سیب 
شهرستان ارومیه پرداخته است. یافته‌های پژوهش حاضر نشان داد که 
همبستگی و معادله رگرسیونی معنی‌داری بین دو متغیر دمای حداقل 
به‌دست‌آمده از مدل ۷۷۲۵۲7 در کل ایستگاه‌ها و مقادیر مشاهده‌شده در 
ایستگاه‌ها وجود دارد. يافته دیگر این پژوهش, قابلیت تصاویر 
ماهوارهای لندست ۷ و ۸ در تشخیص تغبیر مراحل فنولوژیکی درختان 
سیب را تأیید نمود. 

هرچه پیش‌آگاهی سرمازدگی زودتر انجام شود کشاورز بهتر 
می‌تواند برای مقابله با آن آماده شود به همین دلیل در این مقاله 
شبیه‌سازی ۳ روزه دمای حداقل مدنظر قرار گرفت و آرزیابی برای کل 
بطق قمال غرب انران اتسام شف معمول بهدلیل آییکه فیاگی‌های 
سینوپتیکی داثماً در حال تفییر هستنده پیش‌بینی‌ها با افزایش 
مدت‌زمان پیش‌بینی صحت کمتری پیدا می‌کنند ( ۷281627-۷10۳۲0۳3 
8 , .۵ 2۶) در این تحقیق نیز شاخص‌های مورداستفاده برای 
ارزیابی خطای مدل نشان داد پیش‌یابی‌های مدل ابتدا برای 
شبیه‌سازی ۲۴ ساعته و سپس ۴۸ ساعته دقت بیشتری نسبت به 
شبیه‌سازی ۷۲ ساعته دارد ولی با این حال دقت آن با شبیه سازی‌های 
۴ و ۴۸ ساعته اختلاف چندانی ندارد. مولر و همکاران ( ۶۶ ۱۷۲116۲ 
6 ,.2) نیز هرچند کارایی مدل ۷۷2۳ را برای شبیه‌سازی هوای 
۷ روز آینده تأٌیید کردند ولی به این نتیجه رسیدند که همبستگی و 
خطای مدل از پیش‌بینی ۱ روزه به سمت پیش‌بینی ۷ روزه کاهش 
می‌یابد. نصر اصفهانی و همکاران (2019 ,.۵1 6 تصدطدتوظ عدآ) 
نتیجه گیری مشابهی داشتند آن‌ها به این نتیجه رسیدند که صحت 
پیش‌بینی ۲۴ ساعته نسبت به ۴۸ ساعته بیشتر است. در شمال غرب 
ایران که یک منطقه کوهستانی است. شبیه‌سازی جریان هوا در 
نواحی با توپوگرافی پیچیده بسیار مشسکل است و به همین دلیبل 


درمجموع کل ایستگاه‌ها ضریب همبستگی در هر سه شبیه‌سازی در 
حد ۰/۵ و مقدار خطای ۷۸۳۴ و ۷۲5۳ به ترتیب به حدود ۲/۸ و 

محققین دیگری نیز در تحقیقات خود به نتیجه مشابه فوق 
رسیده‌اند به‌عنوان‌مثال سین یو و همکاران (2019 .601 61 جاآ) نیز 
حساسیت دمای هوای نزدیک سطح را به فیزیک سطح زمین نسبت 
به شرایط اولیه و مرزی بیشتر دانسته‌اند. نتیجه تحقیق آن‌ها نشان داد 
که عملکرد مدل در پیش‌بینی دمای نزدیک سطح به پیجیدگی 
توپوگرافی محلی هر ایستگاه بستگی دارد. دوآن و همکاران ( مه 
8 .41 21) نیز نتیجه گرفتند که ارتفاعات و انواع کاربری اراضی 
به همراه تغییرات فصلی به‌شدت بر پیش‌بینی متغیرهای نزدیک به 
سطح زمین تأثیر می‌گذارد. بااین‌حال در تحقیق حاضرء طبق جداول ۴ 
و معنی‌داری مدل رگرسیونی بین متغیر دمای حداقل خروجی مدل 
۳ و متغیر دمای حداقل مشاهداتی درمجموع کل ایستگاه‌ها 
وجود دارد. البته با در نظر گرفتن شاخص‌های بررسی‌شده برای 
ارزیابی پیش‌یابی مدل ۷۷۲۹۳ در هر ایستگاه. ایستگاه‌هایی هم هستند 
۴ الف» ب‌ِ و9 ت). 

بر طبق دومین يافته این تحقیق» شاخص ۱۲۷۲ به‌دست‌آمده از 
تصاویر ماهواره‌ای لندست ۷ و ۸ قادر است پیشرفت و تغییرات در 
مراحل فنولوژیکی سیب را نشان دهد. رشد برگ هم‌زمان با شکفتن 
جوانه در درخت سیب آغاز می‌شود 9 تا مرحله رشد میوه ادامه می‌یابد. 
شاخص نرمال شده تغییرات پوشش گیاهی (۲۷۲() که از انعکاس 
تابش قرمز و ملدون قرمز محاسبه می‌شود. به ویژگی‌های پوشش 
گیاهی شامل شاخص سطح برگ» ظرفیت جذب نور و پتانسیل 
فتوسنتز ارتباط دارد. به‌طورکلی مقدار ۷« با افزایش آشآ یا 
غلظت کلروفیل افزایش می‌پابد (1994 ,عصنته ۷۷ 220 ۷0۵027 ). 
نکته‌ای که در روند مقادیر ۷« در این پژوهش دیده می‌شود این 
بافصل تابستان در ایران مقادیر بالاتری از شاخص ۲۲۷۲ 
به‌دست‌آمده است ولی از سال ۲۰۱۳ تا ۲۰۱۶ بیشترین مقادیر در 
فصل بهار (اواخر مارس تا اواخر ژوئن) بوده است. مطالعات دیگری نیز 
به نتایج مشابهی رسیده‌اند به عنوان مثال نتایج تحقیق خیر و 
همکاران (2019 ,.۸1 2 1611270) پتانسیل ماهواره سنجش‌زدور برای 
تعیین زمان درست رخدادهای فنولوژیکی جوانه با استفاده ازسری‌های 
زمانی ۱۲«۷۲ واسنجی شده در مناطق وسیعی که دور هستند یا 
دسترسی محدودی دارند را نشان داد. ژزنگ و همکاران ( ۶ 70602 
6 ,./0) با ترکیب شاخص ]۱۲2۷ تصاویر ٩۳015‏ و مودیس» 
نتیجه گرفتند که پارامترهای فنولوژیکی به میزان زیادی با داده‌های 


۲ نشریه آب و خاک. جلد ۰.۲۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


زمینی همبستگی دارد. برخلاف نتایج گفته‌شده ون بنویتز آلوارز و 
همکاران (2018 ,۵1 2 ۸۷2۲22 56612 ۷۵۲) به این نتیجه 
رسیدند که شاخص‌های پوشش گیاهی به‌دست‌آمده از 1200926 
+711۷ نمی‌تواند بین مراحل فنولوژیکی متوالی که در آن اندازه برگ 
فقط کمی متمایز است تفاوت قائل شود. البته در تحقیق حاضر سعی 
شده تا با استفاده از تعناد تصاویر و سال‌های مطالعاتی پیشغر دقت 
بررسی مراحل فنولوژی و نتیجه‌گیری در مورد قابلیت تصاویر 
ماهواره‌ای لندست افزایش یابد. 


نتیجه گیری 


این مقاله در پی ارزیابی یک سیستم هشدار سریع سرمازدگی در 
باغات سیب شهرستان ارومیه بر مبنای شبیه‌سازی ۲ روزه دمای 
حداقل و نیز تشخیص مراحل فنولوژی سیب از طریق تصاویر 
اهامای هر تال ها کذشه اس هر ارتیم شام وهای هن 
۳ در شبیه‌سازی دمای حداقل ۷۲ ساعته بررسی شد. مدل برای 
شمال غرب ایران با تفکیک ۴ کیلومتر مطابق با زمان گلدهی سیب 


شاخص‌های ارزیابی. کارایی مدل ۷۷1۲ را برای شبیه‌سازی دمای 
حداقل در منطقه شمال غرب ایران نشان داد. در گام دوم قابلیت 
تصاویر ماهواره‌های لندست ۷ و ۸ در تشخیص مراحل فنولوژی سیب 
مورد بررسی قرار گرفت. با مقایسه ۱۲۲2۷1 به‌دست‌آمده از تصاویر با 
مراحل فنولوژی ثبت‌شده در ایستگاه هواشناسی کشاورزی کهریز 
مشخص شد که تصاویر مذکور می‌تواند مراحل فنولوژی سیب در 
ایستگاه کهریز ارومیه را شناسایی کند. 

این تحقیق کارایی مدل ۷۷۳ برای شبیه‌سازی ۷۲ ساعته دمای 
حداقل و نیز پتانسیل تصاویر لندست ۷ و ۸ در تشسخیص مراحل 
فنولوژی سیب در محدوده موردمطالعه را نشان داد؛ بنابراین با 
به کارگیری خروجی شبیه‌سازی دمای حداقل ۷۲ ساعته مدل ۱۷15۲ و 
تشخیص مراحل فنولوژی از تصاویر لندست, چنانچه در هر باغی با 
توجه به مرحله فنولوژی, دما در ۷۲ ساعت آینده به حد بحرانی برسد 
می‌توان هشدار سرمازدگی را اعلام کرد تا اقدامات مربوط به 
پیشگیری از سرمازدگی توسط کشاورز انجام شود. البته به‌منظور کسب 
نتایج بهتر و دقیق‌تر پیشنهاد می‌شود در تحقیقات آتی از روش 
استفاده‌شده در این پژوهش در باغ‌هایی با مساحت خیلی بزرگ‌تر 


اجرا و برونداد دمای حداقل مدل با مشاهدات ۵۱ ایستگاه سینویتیک 


شمال غرب ايران مقایسه شد. با وجود کوهستانی بودن منطة استفاده شود. 
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۶۲ 19۵6 تما 24 فافهمع0۲ ام0ممه ۷۷۲۲ ۵۶ جمتاقمتله۲۷ .2016 ,۱۳۲ هطاتعظ مه ها مصذاما .۵۷ عملانت۱ 
1001-۰. ۰ (3103 مصتافههع۲۵۲ . مصه .. تعطاهه ۷۷ بجمزعمصده .امد ممطانامد ...عیام مصترماتصمص... متقصصتاهم۲ل بط 
0۰ ۸۳-۲-15 75.۷۷ ۱)105://001.0۲2.10.11 

طططا .عمرمص1 ممتادامعع۷ ۵۵0و ۵۶ طمناها1۳160۲6 16 ,1995 اه علقطونع۱۷ 20 ر.اظ وممااه .۲6 ۲۱21 ر.ظک تممموا۷ 
.481-۰ :(33)2 عصتعمهعو مامصصف 220 0605016866 ۵ 12882661088[ 

0 1 05تقطهته 0عمحصلح عمعصع‌صنهعم عصتادمظ ۵۶ وتومصع۳۵:ظ .2019 ,۴ طقمهم‌صهل ۷22 هه رب۱۷ تصمطهگوظ تودل 
۰( طوناوص ‏ طاذ ‏ صهتوتعظ . ه1) . ,497-512 :(3) 51 . مهتم طمتهمعم . رطامحتومعن . ل2ع1وظظ ,۲۵۵108 
6۰ ۳۳۲0۲۰2019 [.10.22059 8۰ 0://::01.0۲) 

مدمه 107 امه ۷۷۲ ۵ ومتامتل۷2ظ .2018 ,۷۲۸ تصقطحکوظ دا صرح ربظ طمممم‌صه ۷22 ربا تصمطماعظ عوهل 
-163 :(51)1 متام طمتهعوم بطاممتومعن0) آهمتووطاط .صتفحظ بیط طععصعوهم طا عمممععه. مصتاوم0 240 )۲۵۲6628 
۰ ۳ [.10.22059 .۰.001.0۲8 005://0 .(990200ه طفتاعصط ۷160 صهتوع۳ م1) .182 

۶ ط0جمحصمصمطظ )و۲۵0 ومد عوهمتاه مصل2رلممه مه مصاول‌نتاد ,2013 .2 20010صظ مقطعطهنا مه ریک کهتص) 
237-3۰ :(10)19 ۴۵5620 ]مممصم 1۱۵۷۵۱0 لمممتعمک مه امومع ۶ متام[ .ععهز۵۷ظ ۷۵20 ها وهعقطمع) متطمه]18ظ 
(امجتهده9اه طوتاعصط طا۳ صهتوعظ م1) 

۲۲۵۷106 تون عصوزنههاعم2ه ۵۶ صمتاهم2نصدع0) عصتصمهاظ . م4صه . عصمصم‌ممصه] بصمت)ام2تصهع:0 . )عز0باظ . صه . صقاظ 
0۵۵۲ 1۱ ۷151/60) ۳000۲۵1۲ .هدع ۸۵۹://۵2ظ 26 مااه۸۷۵11 .2018-ععص۱۵۱1ظ صو‌ژند۵ ۸26۵۶ ۷۷۵۵۲ 0۶ 001 تهع۷ ل2عاو211) 
2020(۰ 

5۲ ۹688018 (06ص1 ممتاما۷۵۵۵ 01116768686 11260مصصرمظ 0موده-«ممصنهن) ,2013 .لک 0عدهاداب۴۱ 20 ر.لل ۷2۲۵۵ ر.ظ1 همه۴۵ 
135/۰ 2870۳[20. 134 8.10.2 ۳)۸۵۹://01.0۲ .1345-1354 :(105)5 امصتاهژ «حجمممتعخض .عتهاه ممعهعانط جمتامی عطاته] ۳2۵ 
۲12188 0۲۵2۲ عطا ما ممادرد ممتاهامعع۷ مصل:۱۷]۵۵۲]۵ ,۱974 ۱۷۷۰ ع۲۱(6۵:1۵ صه2 .هل 611ظ90 ,الک فمقل ,. ۲۷۷ 6وتام۴ 
309-۰ :351 مم0هع۵۵11ظ م۵11عع(5 مها .۳1 

۲ 10۵ 2)100منام 2 (ولاآه) ۱۱6۱0 اهنوا ۸0116121 2014 ۵ )مهن مصم .۷ طقلطم‌صطوع 1 ره تاتلمطک ر.ظ 2 
فصح 0۱۱0۷ تمانهتتوض ۵۶ لیام تمصمتامصعاص! بعاقل متتقصصتاه عصتفیا معط م1 فععفاه لهعنعمامممطم ۵۶ عمصمتناعع0 
۰ 15739 .6۵://:,:001.0۲8۰10) :352-361 :(2)10 6و۴ 

6 م۵ عصلوت مفمعصف همم 4صقاعوی فصق1 ۵۶ وتوباممم :2013 بک صهتامطوعنا مه ر.ط لملهزامطم 
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۴ نشریه آب و خاک. جلد ۰۳۶ شماره . فروردین - اردیبهشت ۱۴۰۱ 


:(8)37 ظ6210وم اتمه ۵۶ افصساوژ موعتظم .ع0مطامصظط ممتاهم‌لومهاه مه عملم‌ممممنل ( ۷ ۵۱۱۲) «ع1.4 ممتاماهعع۷ 
۰ 7.۸1 ۳1005://001.0۲2.10.589 .4614-4622 

وم .ل ۱10166900 ر.ل 2و1 رم 20 ۷۷ و.ل ملامعت۷]۵0 ,یبا امک ریق تعاطهاو وک عصقطا2 ری مهو .لا تقاط 
۲ 2500۳00۷ ععصه]1۵1160 جوز 500 - ۱۷۲۵۲5 ندرا ,2013 .۳ وم12۷1 20 ریک 9۳2021 ربع2 نا رل ۷۷2۵8۵ .با 24 ,.باظ 
8 مامحصم ۵۲ امصسام[ اقطمتاعصجما. ,ومتصمطر ‏ رعمامممدام ممتاهامعع۷ مصتاماتومص. تم مامتالمم رظن )فمع0ه۵عه 
۰ 2.10 ۳105://01.07 .5997-6016 :(34)16 

6 ۵96۲۷۵0۵ مطا عصلودا واع0مصه عصمصه‌متنوع۲ التطه تمیصته ۵ ممتاجناه۷ظ :2016 هه عه 20210 هه .۴ 21290601 ۲۷ 
۳۵۳۹128 م1) ..561-570 :(4763 ممممزمه اهعلنهمتنم۲۱ مهتصه .رصوی یقتصصیتا) 2تعطفک طا قاقل «عمامصمطم عع8 
۰ 102205911۳19 60۹://001.0۲2۰ 209026)(۰ طوتاعط۴ 

۲ 0۳ 000050۲ 0۲ علو1 مطا مد مفهمت186 .2018 ,۷۲ 06)62م1 0ص2 ر.نا معمافتیط وت مینک رما تعلتمصطامو ر.۷ مووه ۲۷ 
مج احتانهتتوخه ععععع0 عنام )فه1 معط مان مصهاع2ازه ما فصمتاهعاه تعطهنط که قفاوم عصتنمد م) قععع )نبا 2۳00 
.2271011۳061.201 [.۳۲۵05://001.0۲8.10.1016 60-69۰ :248 ۱۷۲۵۱۵0۲0102۲ ۲0۲۵5۲ 

۶ و6عهاه اهمزومامممطن عون :2018 ۱۷۲۰ وع2۷10ععظ-0ع8هه) 20 ریک تقامق-ممصه 22 ری ۷2۲۵2 اظ 102«مصصوظ هم ۷ 
65-1۰ :(36)1 ۳۲۲۸نا] .ومع4ط1 جمتاهاع۷6۵ 20011106-0611۷60و 20 فممتاجموهان 610 تما رما اه مسر تفگ 
۰ -10.4067 .01.018 ,۲0://:۰ 

40 ۵962۲00 تعطهه ۷۷ معط عصونا .2018 ,1 02۲۵20 200 ربکا 200قامن۱ ریت۱2 مهن ر.ل فقطمتلن) ریت ۷1167-110۲0 
۷۰ 1 7/1 عمتعمک صفعهمه صه طا مصافههع۳۲۵۲ ممتاهااض۳۲:۵۵1 ۲۵۲ ۱۷۲۵061 (۷۷۴۳۴) عصتافهد6ع۳0۲ 
47 .2 1005://001.0۲ .304 :(9)8 ع۲عط0ومصص هر 

عصتودا عاممج 0۶ وعفقطم لهمنعم1مصمطام اصمنم۶ن مدتافه۳۵۲6 .2010 ۷۲۰ عهداه)تصمصتعزمو مه .نا تقطام ربظ طفمهم‌مه0 ۱۷۲22 
امناودط طا صعتوع۴ م1 .97-106 :(60۵2 ومممزمو اتمه صا دمتهعومع ۵۶ لقصتاوژ بعاعممه. تهتنامطر. 2101110121 
(ا209020 

220016 ۳۲1۲-و2اع هب1 امه ۵۶ رع4ص1 ممتاماععع۷ 0111670866 ۵۱۴۱۵11260 م1 ,1994 .۴۳ مصلته ۷۷ 220 ر.1ظ ۷۲۵06۲ 
-1]05://001.012.10.1016.0034 .81-91 :49 )صمصمماطظ ۵۶ عصتووعو مامصصفگ .عصمتاه‌اصوعصمی آاوطامممل ‏ مصارونه۷ 
94(90061-۰) 42517 

عط صا صد رن ممتاهاامهعم 10۲ فافهم1016 ۷/۴۲ 0۴ ۷۵۳۱6۵20۵08 .2014 .۲ موی لحم ری ت22 رم م2 
رامجمتوطاه طوتاعص از صعتوع۴ 12) .3-10 :(38)87-86 ,۱1۷۵7 :2009 ۳۵۱۵-۱۷۲۵۷ 6۳100 

9۵200-1۵۵ طعناط مصلونا ممتاهعاه۱ «ومامعمطط مومت ,2016 .۱ عصع2 4ص2 رب عصقط2 ر.ظ ۷۷۲ .۷ 2۳698 
«(16)12 ,5605015 ۳۲0۰ ۱/9 230 9۱۳15 11۳01 و۱۱ 12 ۵00 
9 ۵/0102 2099 

۶ م۳609 صرح محتناهعه وم ا06ممهر ۷۷ ۵۶ همتاهم۵۳1/1 ۷ .2013 .2 1تدگل هه ,۷ ۱۷۲۱۲۹۵۵۵0 زمهدان ربا ت0وز701 
8 .187-194 :(2802 اهوم لهءنامدتوهع0 صویا صا اقهمع۲۵۲ ممتامالمههين ۵ ممتمتلهبه ده ومصصفطعو فنامته۷ 
امحتاوطاه طفتامدظ طزه مقتو۳۵ 
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